function [procX procT] = signalPreparation(signal,eps)
%**************************************************************************
%vectorization of the image

%input:   path to the image file
%output:  Matrix X which is NxT type where T=PxQ where P and Q are
%dimensions of image in pixels and N=3

%**************************************************************************
dim=size(signal);
T=dim(2);
procXCal=zeros(dim(1),T);
count=0;
for k=1:T
    if norm(signal(:,k))>=eps
        count=count+1;
        procXCal(:,count)=signal(:,k)/(norm(signal(:,k)));
    end
end
procT=count;
procX=procXCal(:,1:procT);
clear procXCal;
clear signal;